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NATIOML ADVISORY COMMITTEE FOR AERONAUTICS. 

TEOHRIOAL MEMORAtE)UM NO. 325. 

COMPUTATION OF CANTILEVER AIRPLANE T/INGS-* 
By K. Thalau. 

The so-called cantilever wings, often preferred in recent 
times, introduce, in conjunction with braced girder construc- 
tion, high-grade statically indeterminate systems into the com- 
putation. We will first consider a cantilever wing with two 
spars and then a similar v/ing with three spars in an approxi- 
mate computation. 

The purpose of this treatise is, first of all, the deter- 
mination of the effect of variously loaded spars on one an- 
other, since the neglect of this effect would present an eco- 
nomically very unfavorable computation method (See also the 
article by L. Ballenstedt, Technische Berichte, Vol* III, No, 4). 

The system of spars and cross-bars alone (whether solid or 
built-up) does not matter at first, the original assumption 
being that, the spars are rigidly braced by the cross-bars. 

Two-Spar T/ing 

Our system can be regarded as a one-sided Vierendeel gird- 
er, on which the principal vertical components of the air 
forces act perpendicularly to the plane passing through the axes 
of the spars at the junction points (Fig. l) . 

* From "Zeitschrift fur Flugtechnik und Motorluf tschiffahrt," 
May 26, 1924, pp. 103-109. 
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The magnitude and direction of the forces F and F' change 
according to the manner of loading, hut are to he regarded as 
constant for a given load. Hence the load diagram of the air 
forces in this direction, contrary to their actual distribution 
over the width of the wing, is assumed to be rectangular. This . 
assumption corresponds to the method of computation often emr- 
ployed. We are independent of variations in the dynamic pres- 
sure distribution,, when we make the computation for the general 
case of a spar loaded with unit forces at its junction, points. 
This produces certain tt forces or moments, which we will desig- 
nate by A". An m-fold load (m x 1 = F) would accordingly pro- 
duce an mA' force or moment. An n-fold loading of the rear 
spar (n X 1 =,--F') gives a corresponding symmetrical result of 
the form A = mA' + nA'. 

We will now give the computation for a system with five 
fields of constant width X and constant height . 2h. We get our 
statically determinate main system in the form of two fixed gird- 
ers by cutting through the cross-bars (ribs) and obtain, accord- 
ing to Vierendeel (Fig. .2), by combining the unknown lateral 
forces,, the following new unknown quantities! 





TT, + TT, + TT^ + TT. 
13 3 4 



+ + TT3 




= X; 




3 



X3 



(1) 
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In the case 2X = 0, the front spar is acted on only hy 
the external loads "1" . . 

In the case = - 1,. all other X values (in the compar- 
ative group l) equal 0 and we obtain TT^ = - 1^ TTg = + i. 

Oh following up the different cases of X = - ly we accord- 
ingly find that always two of the tt forces are effective, which 
together give simple hending and torsion moment surfaces. An 
exception is formed by the case Xg = - 1, through which only 
TrB= - 1. 

As follows from the subsequent figures, bending and torsional 
moments for the spars and bending moments for the cross-bars 
(ribs) are the most important in the elasticity equations. 

Disregarding, therefore, . the perpendicular and lateral 
forces^ the general expression for the bending deflections is 

5., = f ^i ^ r Mj Mk dx r Tj T], dx 

^s •'•s % •'•V 

in which signify: 

Mj[,M]5., -bending moments for the spars; 

" " » " ribs;: 

Ti,Tj5., torsion " " " spars; 

Ejjj elasticity modulus of a spar (same for all spars); 
Eg, " " » " rib ( " « " ribs ); 

inertia moment of a spar ( « n ti spars); 
Ig, " " II rib ( " " " ribs ); 

Gh, gliding modulus of a spar (same for all spars); 
Ip', polar moment of inertia for a spar with reference 
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to the deviation of its cross-section from the circular shape 
(same for all spars). 

The relatively small angle of -torsion of two cross-sections 
1 cm apart (Hiitte, Edition 23, Vol. I, p. 570) is, according to 
Grashof, 

G = £ — hm % = . . £ ~ 1.2 

* ^xH lyH ^pH« % 

If, starting from constructions at hand, we express lyjj 
"by. IxH> "then, with lyjj = and Ip = + ly, we obtain 

= ^xH + = ^xH ( ^ + a / 

We then obtain 

T ^ ^xH ^xH _ ^ T 

a (1 + ^ A (, I „ 1.2 a ( 1 H- •=!^ > 
V ay xH V a/ 

^pH' " a^+^1 ^.xH = ^ ^xH 

in which b is mostly- a true fraction. 
Moreover^ if we put 

^s = ° 

Gh = d Eh; 

Eg = e Eh; • 

so that^ after substituting these values in equation (2)^ we find 
for Siic only the elasticity modulus and the inertia moment of 
the spar 
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Mj dy. 

Eh Ih 



+ / 



Ml 1% dx 
e Eh c Ih 



+ / 



Tj, 

d Eh "b Ih 



With 



e c 



c' and ^pg- = c" ■becomes 



■^i dx 



+ 9" / 



dx 

Eg Ih 



(3) 



The elasticity equations of our system now read 



^0 2 



= + Xg 633 + X3 .633 + X4 634 + Xg 5; 



^©5 = Xi 651 + Xg 652 + X3 653 



+ X4 5b 4 + Xs- S55 



J 



(4) 



By Multiplying 'botia, sides of the elasticity equations by 
Eh X Ih» "fcT^ese denominators are entirely eliminated* 

We will now take up the representation of the moment areas 
and 'their evaluation for the different cases. 

Case 2X = 0. - Only the external forces "1" act on the 
front spar (Fig. 3) . 

Tjg = 1 X 

•n^ = 1 (2\ + \) = 3X 
Tia = 1 (3X.+ 2\+ \) = 6X 
Tig= 1 (4\+ 3K+ 2X+ \) = 10 \ 
Ti^= l {5X + 4\ + 3\ + 2X + \)= 15a. 

Oase X^ = -1. - We have tt,^ = - 1; tt^ = + l (Fig. 4) . Mo- 
ments which generate pressure in the upper fibers of the girder 
are positive. The bending moment is the same for all ribs and 
equals 1 X h. At the same time, there is a constant torsional 



P = ^ x^ 
P, = 1^ X^ 



(5) 
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moment for the spar s> whereby it should be noted that the latter, 
with the exception of the case = - 1, only acts within one 
fieldi with the case = - 1, only for spar sections in field 
2; with case Kg = - 1,- only in field 3, etc 

The torsional moment areas , on account of their simplicity, 
were not plotted and are easy to follow.' 

Case Xg = - i ; TTg = - 1; TTg = + 1 (Fig. 5) * 

Gases X-a = - 1 and Xa = - 1 give analogous moment areas. 

Case Xs = - 1 .- This gives TTg = - 1 (Fig. 6). Here the 
influence of the torsional moment extends over the whole spar 
length with the magnitude 1 x h. 

Regarding the evaluation of the simple moment areas, see 
Mtliller-Breslau, Voli II, Chapter I, as also the handy formulas 

of Demel«* 

* Richard Schadek von Degonburg and Karl Demel "Hilf smittel, " etc., 
Berlin, 1915, p. 7 - Published by Wilhclm Ernst and Sohn. 
If, e-*gi, we wish to find the value of 

= / ^3 dx ^ J ^1 ^x ^„ J Ti T3 dx 
Eh Ig 3h Ih Ejj Ih 

the moment areas for the cases = - 1 and Xs = - 1 are to be 
combined. (Figs. 4-5) i The integrals are to be extended only over 
the portions of the bars which are simultaneously subjected to mo- 
ments'. Hence ^ 

/ li. Ma dx = / (+^) ( + X)dx + / (-X) (-^) dx + 



^ Field 1, 3^ Field 1,. 



behind ^par. ^'^'^'^ spar- 

.+ /^ Xx dx +/ (->) (-X) dx 

: Field 2,, - y y Field 2,. 



behind spar. front of spar. 

2 



• / Ml -Mg dx = 2 )? + 2 p = 3 \^ 



c' / Ml Ms dx = 2 c / (+x) (-x) dx = - I c' h^ 

1. ?Cr OSS-bar 2 



(Continued bottom Page 7) 
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Demel' s formalas give: 
For the load magnitudes - . 



£ 




=. / 


Mo d X = 


-101 « 
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I §0 3 




23 t 
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z* 


£ 




= + 


100 X® 






Z5 



For the form magnitudes - 
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I 


5x1 


8 
3 
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3 °' 


If 
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I 


5l3 


= 3 


t 




3 °' 


h^ 
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613' 
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c" 


h^ 


X 



(Continuation of footnote from Page 6) 

c" / T,^ 1^ dx here equals zero. Hence 



The work of integration is rendered unnecessary by employing 
the above formulas. The moment areas are, due to the arrangement 
of the unknowns, almost all rectangular or triangular surfaces, 
the torsional moment areas being also rectangular surfaces. Trap- 
ezoidal surfaces are divided into triangles by means of diagonals. 
The signs must be given attention. 
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Ej^Ssi = E I 6^2 

E I 633 = :C + I c« 1?+ So" li X . 

E I £te3 = 5 >f - 2 c' ^ 

E I 634 = 5 ;f 

E I fifes = - - 2 c» li \ 

E I = E 1 5^3; SI ^3 = E I &3 

EI 2^3 - ^ ^ 3 °' + 2 c" if X 

EI 6&4=7 X^ - 

E I 635 = - \^ - 2 c" f \ 

EI = E I 5i4; E I 64-2 = E I E I 643 = E I 6^4 

E I = ^ f + I c' 1? + 2 c" ^ X 

EI 645 = - ^ X^ - I c' - 2 c" li X 

EI 651 = E I E I = E I 635; EI 653 = S I &3s; 

EI 654 = E I 645 

E I 655 = ^ ^ + I c' + 10 c" f X 



We designate the load magnitudes "by to. Zs 
the elasticity equations in the form of a matrixi 



and write 
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X3 


• 








5aa 


5.3 




5x5 


Zs 




5a s 


53 g 


5a4 


535 


Z3 


531 


63 a 


633 




63s 


z* 


5*1 


6« 


643 




5*s 


Z5 


56 1 


6s 2 


553 


654 


5s s 



If we iATish to follow the effect of the loads on the differ- 
ent unknowns > we thus obtain, for sLhy gll/en unknown, the form 

Xi = Pii Z^ + ^3 + Z3 + 4- 

For this purpose the corresponding matrix reads: 





Zx 


Zs 


Z3 


z* 


Zs 




3xx 




PX3 


PX4 


Pxs 


X3 


Pax 


Pas 




Pa 4 


Pss 


X3 


fei 


^3 


P33 


^34 


P35 . 




P41 


P43 


343 


P44 


P4S 


Xs 


Psx 


Pss 


P53 


P5 4 


Pbs 



Since = 5-|^^, then '^^■^ must equal i.e., the 

p matrix, like the 6 matrix, is symmetrical to the main diag- 
onals (underlined). The P values may be found simply by the 
methods of Muller-Breslau, Vol. II, Chapter I, p. 178. 

As a summary, we will give the solution diagram for this 
special case of five elasticity equations with five unknowns lii 
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in each equation. From the general case of the niiie-membered 
elasticity equations, v/e obtain the following: 



r 


a^ 




Cr 




©r 




C;j* 1 


b-j,i 


a^-i 




Mr 






^r 


1 










5ti 


5,2 


5l3 


614 


5i5 


^1 
^11 






^1 




3 








5tei 






5a4 






""l 
612 


P2 


^8 






3 






531 


63a 


533 


534 


•53s 






1 

6x3 










4 




6« 


5« 


5« 


5*4 


545 








1 










5 




fibs 


653 


5s* 


55 s 








1 


1 
©IS 











Taken in order, the coefficients are: 



Oil °11 °ll 11 

ei3 = 5^3 + cLis 1^1, in which 
<3-is = 5gi. 







+ 


533 














+ 


fis4 






^13 






tl 


+ 


5.5 



^13 = 533 + C„ + ti2> in Which 

Cl3 = 53^ 

dja = 633 + 631 Hi 

tr 
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©I* = 



3 4 



Cl3 <?! + ^13 ^3+5 
C„ + d.13 ^3 + Ass 



'13 



+ «Pi + Ci4 ^3 ^14 1^3* in Which 



1=14 = 



'14 



M-4 



T3l4 ^1 + % + 1?3 + 643 



e 



14 



'15 



655 + 553. + his (?a + Cjs + «ii5 M.4, in which 
T3is = S33 + 6s 1 

Cl5 = 653 + ^Bl ^ + ^5 

d-lB = 654+ 5fel + 1^15 ^2 + ClS 1^3 

Beginning with the last ^ values, we proceed to make out 
the P tables. 



A Pl5 




^ Pas 


+ 




P35 


+ 


^1 ^^B 


Ps5 




Ms Pas 


+ 




P45 


+ 




P35 




M!3 ^43 


+ 


^3 








^45 




^4 ^55 






























®1B 
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^ Pa 4 


+ 




P34 


+ 


^X 


P44 


+ tx 


P54 




=■ 


P34 


+ 


^2' 


P44 


+ 


^3 


P54 






^34 


rr 


Ms P44 


+ 


V3 


•P. 4 
















M'4 Ps4 


+ 


1 

^4 












^S4 




P45 . 


















A Pi3 




II A 

H-i Pas 


•4- 


X 


r^3 3 




^X 


P43 


+ -vlr 


P53 


Pp., 

~2 3 




S '33 


+ 


V 

3 


B 

1^43 


+ 


•fa 


Ps3 






Psa 


— 


Ms P43 


+ 




Ps3 


+ 


1 






^43 




P34 


















PS3 


= 


P3S 






















^X I^SS 


+ 


^X 


D 

Psa 


+ 




P4a 


+ ti • 


Psa 


Pa 2 


_ 


Mis Paa 


+ 


^a 


P4a 


+ 


^a 


Psa 


X3 




"^3 a 





Po, 
^^3 3 


















"42 





•^3 4 


















•^5 3 





•^3 S 






















Mx Psx 


+ 


^x 


P31 


+ 




P4X 


U. -(If 


0 

Psx 


Pax 




PX3 


















Pax 




Px.3 


















p4X 




PX4 


















Psx 




Pxs 



















The arrows indicate the order of the solutions. The unknowns 
can now he computed from 8, whereby no further difficulties inter- 
fere with following, in the manner indicated at the beginning. 
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the effect of the variously loaded spars on each other. 

Three-Spar ¥ing 

From what has been said, the investigation of the three- 
spar system is now simple (Fig. 7). As before 



X. = 


^1 


























Xx 




Xs = 




+ 
























Xa - 


Xx 


Xa = 


■^1 


+ 




+ 


^3 


















X3 - 


Xs 


x* = 




+ 




+ 




+ 












TT4 




X4 — 


X3 


Xs = 




+ 




+ 




+ 




+ 












X5 - 


x* 


Xe = 


^6 


























Xs 




X, = 


^6 


+ 


TT, 






















X7 - 


Xs 


Xe = 


ITS 


+ 


TT, 


+ 














\ 






Xs - 


X7 


Xg = 


TTs 


+ 


TT^ 


+ 




+ 












TT9 




X9 - 


Xs 


X.o = 




+ 




+ 






^9 


+ 


^10 . 








Xxo- 


Xs 



The TT forces are at first assumed to be produced by loading 
the. front spar alone with the junction-point loads "1". The' ef- 
fect of the forces' F = m x 1 then equal mA' . The forces 
F' = n X 1 in the rear-spar junction points again produce a sym- 
metrical or bisymmetrical result. 

'If we now cause the loads "1" to act in like manner on the 
middle spar, then, because of the symmetry of the system, the un- 
knowns, thereby produced, are ir^ = TTg, TTg = tt^, etc. This will 
be manifested by the production of 10 correspondingly symmetrical 
6 values. 
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Since Xj^ = Xg, X^ = X,, etc., only one of the two equation 
groups is necessary for making this part of the computation, 
whereby this part is immediately? carried hack to the case of the 
five equations with the five unknown quantities. If this last 
loading case, with the Junction-point loads r x 1 = F" , has an 
effect r x a"« on our system, then the total effect, when the 
forces F, F« and F" act simu:.taneously and respectively on the 
front, middle and rear spars, will be 

A = m A' + n A" + r A'" 
We now have at our disposal 10 elasticity equations with 10 
unknowns, as follows: 

Sqi = Xi5^^ + Xs 6.^3 + Xa + X4 + Xs 615+ 

+ Xe + X, 5,, + X3 + Xg, + X,o 5„„ 

+ Xe + X, 63, + Xs + Xg 639 + Xio d^^^ 

m 

Xe 610s + X, 5io7 + Xg 5^0 8 + Xg 5^09 + X^Q Sjolo* 
The general structure of the form magnitudes fij^j^ is the 
same as already given on a previous page. As Fig. 8 shows, the 
load magnitudes 6q^ to . equal zero, while to 5^xq ful- 

ly agree with the values of 5^^ to S^g already given. 

The combination of the Mq surfaces with the - Mg sur- 
faces make, as already mentioned, 6^ ^ to 6^5 equal to ziero 
(Figs. 9-11). 
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¥e then have 

■r, P 101— -jS . 

E I bof ==— 23 ; 

E I . 5 o 8 == - - ^ '> ^ 

E I 6o9 — ~ ^ > 
E 1 . 5a 10 = + 100 ?? ; 
The form magnitudes read: 

E I 6ii'^ E I Sj- 

E I 6is = - I ^ + <5" 1^ ^ 



E I 5i7 = - J 



2 



E I 6.xe = - I 



EI 5i9 = - I 



E I 6110= i?- \^ - p'" ^ 

E I 52 1 E 1 \s 

3 
2 



E I - 5^6 = - I t {= 61,) 



E I Sgv = - I + c" la X 
E I fee = - I 
E I fe g = - I 



E I 52 10= + ^ X - c" h \ 
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I 


631 


E 


I 635 


E 


I 


Sas = 


— 


1 ^ ( 


E 


I 


63, = 




f ^( 


E 


I 


^3 8 = 




mi 


E 


I 


639 = 






E 


I 


^.3 10 = 


4' 





'18) 



E 


I 




E 


I 645 








E 


I 


64a = 






( = 






E 


I 


S47 = 






( = 


639) 




E 


I 


648 " 






( = 


§3 93 




E 


I 


5-4,9 = 






+ 


c" 


If X 


E 


I 


5*10^ 


+ 






c" 


^ X 



E I 651 E I 653 

E I = + ^ \3 - c" 1! X (= 6,, J 

E I 5s7 = + ^ ^ - c" if X (= 62x0) 

E I = + ^ - c" h= \ (= 6310) 

E I 659 = + c'.' -^^ .(= 5*10.) 

E I 6510= - + 5 c" li X 
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§61 = 


E 


.1 5l6 


E I 5,, = 


E 


I 6,^ 


§63 = 






6?3 = 




^7 


^63 = 






673 = 




^7 


5s* = 






674 = 




647 


.5e 5 =" 




656 


67B = 




657 


e 

°ee = 




611 ■ 


07 8 = 




6s 7 


6s 7 = 






V,/y 7 — 






§6 8 — 




0*1 0 






&^ a 
"-733 


5eQ = 




"14 


'-'7 9 




Mff 4 


10 




"IS 


* 10 






Ssi = 


E 


I 5l8 


E I = 


E 


I 6i9 






^8 


6b3 = 






Sea = 




^8 


^3 = 




^9 


Se* = 




648 






649 


Ses = 




S58 


^5 = 




^9 


c — 
°86 - 




0S8 


. ^6 = 




639 


O f 




5k» e 
7 B 






^9 


8 




"3 3 






"89 






fii34 






44 


A — 

«8lO= 






^10= 




645 




E 














5l0 2 = 


6210 










6103 = 


§3 10 










610 4 = 


6*10 
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6, „ = = E 1 


10 


5. . „ = 




0l0 7 = 


0? 10 


Sio 8 - 




^10 9 = 


^9 10 


^10 10~ 


55 5 



If we again use a matrix for summarizing the ten elastieity 
equations, we then oTatain, on designating the load magnitudes 
by Zito Zio, 









Xg 




X5 




X7 


Xs 


X9 


Xio 




0 




5i3 


5^3 




5is 


5ie 


5x7 


5x8 


5x9 


5x10 


^) 


0 






^3 




^5 




§37 


^8 


^9 


^10 


^) 


0 










635 


§36 


537 


638 


639 


5310 




0 










54,5 




5*7 


548 


5*9 


5410 




0 












55,6 


557 


55 8 


559 


5s 10 


^) 


Ze 














5s7 


5s8 


569 


56 10 




Zv 














577 


67.S 


579 


5710 


^) 


Zs 
















§88 


5.8.9 




^) 


Zg 




















59 10 


10) 


Zio 




















5x0 10 



Herein it is now worth noting that (as the computation of 
the Sj^]^' values demonstrated) there is a symmetry in so far -as 

5ii = Sgg; 632= etc* 

5is = 667; • 5i3= 588 etc. \ 
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Hereby it is also possible for us tb deduce the t^n^membered 
elasticity equations to two groups of five equatiblife With 5 un- 
kncras. We first add equations 1 and 6, 2 and 7, 3 and 8, 
4 and 9, 5 and 10, and thjgn subtract them from one another » 
The resulting members are then designated briefly as follows:- 





+ 


Xe = 


Si 


^.1 - 


Xe - 


Di 




+ 




Sg 


Xg - 




Da 


X5" 


+ 




S5 


x/- 


Xio~ 


D5 


5xx 


+ 


= 


61 1+ 


5x1- 


5x6 = 


6x1- 


Sxs 


+ 


Sx7 = 


6x8+ 


• 5,3 - 


6x7 = 


5l3~ 


Sxs 


+ 


5x30 = 


5x5+ 


5x5- 


5x10= 


5x5- 



and both groups of the elasticity equations are obtained in the 
form 





Sx 








S5 




5ii+ 


5x2 + 


5x3 + 


6x*+ 


5x5+ 


^ 




633+_ 


633+ 


5te*+ 


635+ 
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The solution of these eqtiations has already been discussed.- 
■A-ft.er obtaining the values S or D, V7c quicKLy find the un- 
knowns X and from them,, in turn, . the tr forces at the inter- 
section points of the ribs. The latter computations make the 
least work., 

We will now consider briefly the loading case of the middle 
spar. Here,... as already stated, there are ten symmetrical load 
magnitudes (Fig, 14). ¥e will also choose a symmetrical direc- 
tion of the TT forces (Fig. 15) corresponding to the bending 
of the system. 

Hereby, in the cases Xg to X^o = - 1, "the signs of the 
Ms "to Mio surfaces alternate, so that the combinations of the 
Mq surfaces with the Mg to Miq surfaces all have the same 
sign,, negative in this case.- The first five equations are all 
correspondingly to be provided also with a Z " member. 

In conclusion, we note that the computation difficulties 
lie almost exclusively in the solution of the many-membered eqxia- 
tions.. . From this side also, attention is called to. the fact that 
only work with the calculating machine offers, in many places, 
promise of good results.. 

The usTial utilization of determinants is inexpedient in the 
present instance, .where the number of equations and unknown 
quantities exceeds four (Compare Muller-Breslau, Vol. II,. Chap- 
ter I). . 

After the decisive forces have been determined, .they natur- 
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ally give the bending moments along v/ith the torsional moments,' 
which now, however, extend their Jeff ect over the whole length 
of the spars. For strong torsional effects, such as may he pro- 
duced at the base- of long spaxa, the stresses produced in^ the 
longitudinal direction of the sjjars must be considered (See 
Goetzke, Z.d^V.d.I., 1909, p,93|) . 

This treatise is intended, first of all, to combine, in 
the example discussed, the most igssential theoretical consider- 
ations with the practical viewpoints so decisive in airplane 
construction. It is, of course, possible to employ more refined 
methods of calculation, such as, e.g., the introduction of the 
torsional moments and the transverse forces in the girder plane 
at the intersection points of the ribs, as further independent 
unknowns, but the work of computation would then be entirely 
disproportionate to the practical advantages gained. 

A numerical evaluation of the above may be given in a future 
article. 

Translation by Dwight M. Miner, 
National Advisory Committee 
for Aeronautics. 
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Figs.l&2 
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Fig.3&4 




Fig, 4. 
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Figs . 5&6 




Fig, 6 
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Figs.7&8 





Fig. 8, 
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Figs. 12 & 13. 
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Fig. 15. 
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